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We study the thermoelectric properties of three-dimensional topological insulators 
in magnetic fields with many holes (or pores) in the bulk. We find that at high 
density of these holes in the transport direction the thermoelectric figure of merit, 
ZT, can be large due to the contribution of the topologically protected conducting 
surfaces and the suppressed phonon thermal conductivity. By applying an external 
magnetic field a subgap can be induced in the surface states spectrum. We show that 
the thermoelectric efficiency can be controlled by this tunable subgap leading to the 
values of ZT much greater than 1. Such high values of ZT for reasonable system 
parameters and its tunability by magnetic field make this system a strong candidate 
for applications in heat management of nanodevices, especially at low temperatures. 



Introduction. Topological insulators (TIs) recently attracted a lot of attention as po- 
tential candidates for spintronics applications.'^ This interest is based on the remarkable 
properties of TIs,'^ namely, the fact that their surfaces (for three-dimensional TIs) or edges 
(for two-dimensional TIs) possess topologically protected conducting states. In this paper 
we argue that these unique properties of TIs can also be employed to make very efficient 
thermoelectrics.'^ 

Good thermoelectric materialJ^^'^ are characterized by a low thermal conductivity and 
high thermopower and electric conductivity. The large thermopower requires steep depen- 
dence of the electronic density of states on energy, which can be achieved by having the 
chemical potential close to the bottom of a band. The relatively high conductivity demands 
the gap to be low. These combined requirements point to semiconductors with heavy ele- 
ments as the best candidates for the thermoelectric materials. Some of these semiconductors 
were recently rediscovered as TIs. We show that by making holes in these TI materials (see 
Fig. [T| and inducing subgap in their surface states by means of applied magnetic fields 
or proximity of ferromagnets one can significantly enhance their thermoelectric efficiency. 
The high density of holes in the direction of transport has two positive effects: 1) to trap 
phonons and thus reduce the thermal conductivity; 2) to increase the surface to bulk ratio 
and therefore effectively enhance the electric conductivity of the sample. The effect of the 
applied magnetic field is to increase the thermopower. 

Thermoelectric figure of merit. The thermoelectric efficiency is defined as a ratio of 
an electric power generated to the total heat fiux through the sample. Therefore, it is 
characterized by the dimensionless number 

ZT = — T, (1) 

Kg + Kpfi 

where a is electric conductivity, S is Seebeck coefficient (or thermopower), T is temperature, 
and K, is thermal conductivity. The latter, in general, has contributions both from electrons 
(Kg) and from phonons (/tp/i). 

In the absence of magnetic field or magnetic impurities, the TIs have ungaped surface 
states. These propagating states are confined to the close proximity of the surfaces. Their 
existence is protected by the topology of the bulk band structure. These states have cone-like 
2D Dirac spectrum, E = ±vh\k\, where v is the constant Dirac electron velocity. Application 
of a magnetic field, doping with magnetic impurities, or hybridization of these states due 
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FIG. 1. (Color online) A sample of the topological insulator with many holes in the direction of 
transport. 

to close proximity of two surfaces can induced a Dirac subgap Ajl^ZHll -which in the case of 
magnetic field B is A oc B. Then the surface spectrum takes the form E = ±\/v^Wk^~+~K^ , 
see Fig. |2} Generally this subgap is much smaller than the bulk gap, A ^ Aq. For three- 
dimensional TIs Aq ranges from 0.15 eV to 0.3 eV for such as Bi2Te3 and Bi2Se3. 

To study the thermal and electric transport in these materials we employ the linear re- 
sponse theory.'251 The electric (j^) and thermal (j'^) currents are given by linear combinations 
of the temperature and chemical potential gradients: 

f/e = LoVfi + L,{VT)/T, (2) 
f = -L,Vfi-L2{VT)/T, (3) 

where e is the electron charge. Using Onsager relations, one can find from these equations the 
electrical conductivity a = e^Lo, Seebeck coefficient S = —Li/{eTLo), and electron thermal 
conductivity He = — L\) / {TLq). The figure of merit, ZT, can then be represented in 
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terms of these linear coefficients as 



r , ^rr.. r. ' (4) 



Lo{L2 + KpHT)-Ll 

In Eq. ^ it is assumed that the transport coefficients have bulk and surface contributions 
Ln = n + whcrc V = {A — J2n ^-^n)/ J2n '^'^^in is the factor related to surface/bulk 

ratio (porosity) and has dimension of length. The holes in the sample does not have to form 
a periodic structure and can be placed randomly. The same idea applies to porous TI 
samples, see Fig. [3] (a). The parameter T> characterizes the average distance between the 
pores (holes) of the average radius R, see Fig. [3] (b). Here Kph is the phonon contribution to 
the thermal conductivity in the bulk (phonon contribution to the thermal conductivity for 
the surface TI states is much smaller than that in the bulk). 

To estimate the surface contribution of TI to the transport coefficients L„ we assume the 
bands to be Dirac-like with a subgap A and use Boltzmann equation in the relaxation time 
approximation, 



— oo 



Here the sum is over the upper and lower bands, i = ±1, and f'{E) = df /dE with / = 
l/(e*^'^~'^-'/*-^^^^ + 1) being the Fermi distribution function. Then taking relaxation time r to 
be independent of energy, we find 



' ax- 



X — /i)" ^ (—X — /x)" 



cosh^^ cosh2 2^ 



(6) 



2/i2 X 
where h is Planck constant, A = A/{kBT), and Ji = ^/{ksT). 

When the chemical potential /x is far below the bottom of the conduction band, (Aq — 
fi)/{kBT) ^ 1, the contribution from the bulk to the transport coefficients is exponentially 
suppressed, „ oc e^^^'^^'-''^^^''^'^\ and can be neglected. The same argument works also for 
the valence band contribution when /x is far from the valence band edge. Thus, the thermo- 
electric transport is dominated by the surface states and the only sensible bulk contribution 
is to the phonon thermal conductivity Kph. Then the figure of merit becomes 

— ^ 

At small V, the contribution to ZT mostly comes from 2D surface states and in this limit 
ZT is given by ZT2D = -^^ i/(-^s,o-^s,2 — -^^i), which is shown as a function of the chemical 




FIG. 2. (Color online) A sketch of TI's band structure with the subgap A in the surface band. 
The bulk gap Aq is considered to be much larger than A. 

potential in the inset of Fig. |4] (a) for A/(A;bT) = 3 . The color plot of ZT2D as a function 
of the induced subgap A/{kBT) and n/{kBT) is shown in Fig. |4] (a). The maximum ZT 
achievable by tuning the chemical potential for a fixed induced gap A/lksT) is shown by 
a golden line in Fig. |4j The very high values of ZT in reality are not reachable since any 
small contribution from the phonon thermal conductivity reduces ZT. 

For the holey sample, the phononic contribution to the thermal conductivity can be 
characterized by the dimensionless parameter Kp^ = 2nphDh^ / [rk^^T'^). Here the phonon 
thermal conductivity can be estimated to be Kph ~ lWm"^K-i (as for BisTeaP and the 
average distance between the holes reaching P ~ 10 nm. Taking the relaxation time r ~ 
10^^^ s, at room temperatures we estimate Kph ~ 1 for these rather conservative values of 
parameters. In Fig. |4](b), estimated ZT is shown as a function of A and for the phononic 
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FIG. 3. (Color online) (a) A schematic view of a TI sample's part with pores propagating in 
the direction of transport, (b) Top view of the TI sample with the average distance between the 
holes/pores given by P. 

bulk contribution characterized by the dimensionless parameter Kph = 3. Although reduced 
considerably from its pure 2D value, ZT remains substantially larger than any value thus 
far achieved in these materials and can be tuned significantly by applied magnetic fields. 

Summary. We have studied the thermoelectric properties of three-dimensional TIs in 
magnetic fields with high density of holes in the sample. By applying an external magnetic 
field a subgap can be induced in the surface states spectrum. We show that the thermoelec- 
tric efficiency can be controlled by this tunable subgap. We have find that the thermoelectric 
figure of merit ZT in these materials can be much greater than 1. This is due to the high 
contribution of the topologically protected conducting surfaces and the suppressed phonon 
thermal conductivity. High ZT values for reasonable system parameters and its tunability by 
magnetic field or magnetic impurities make this system a strong candidate for applications 
in heat management of nanodevices. 
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for insightful discussions. This work was supported by NSF under Grant No. DMR-1105512, 
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1678). 
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FIG. 4. (Color online) Thermoelectric figure of merit, ZT, as a function of surface subgap A and 
chemical potential fi. Light solid curve shows the maximum of ZT. (a) ZT for the infinitely high 
density of holes. The inset shows ZT{fi) for the case of A/{kBT) = 3. (b) ZT for a finite density of 
holes characterized by the dimensionless parameter Kp^ = 3. Because of the phononic contribution 
to the thermal conductivity ZT is reduced compared to the case (a). 
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